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Gas-liquid criticality in the ultrasoft restricted primitive model (URPM) of polyelectrolytes is studied using the 
collective variables-based theory. For the model, an effective Hamiltonian is derived and explicit expressions 
for all the coefficients are found in a one-loop approximation. Based on this Hamiltonian, the phase and critical 
behaviour is analysed. Our results provide evidence that the nature of the gas-liquid criticality in the URPM is 
the same as in the restricted primitive model that includes a hard core. 
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Recently, a new model of polyelectrolyte solutions, the so-called ultrasoft restricted primitive model 
(URPM), has been introduced in Glia . The URPM is a globally electroneutral equimolar mixture of positive 
and negative spatially extended charge distributions. Contrary to the familiar restricted primitive model 
(RPM) O, no hard cores are involved in the URPM. Nevertheless, the URPM admits a thermodynamic limit 
since it is H-stable in the sence of Fisher and Ruelle HID. The phase diagram of the model was studied 
in Giiiia using Monte Carlo (MC) and molecular dynamic simulations. The analysis of simulation data 
for relatively small system sizes revealed a nontrivial topology of the phase diagram: a region of low 
temperature gas-liquid phase coexistence terminating at a critical point and a sharp conductor-insulator 
transition above this point (ll|2). The simulation data also hinted at a tricritical nature of the gas-hquid 
critical point. Such a phase behaviour qualitatively differs from that found in the RPM. More recently 
(H, the gas-liquid phase transition and the critical behavior of the URPM have been studied using finite 
size scaling MC simulations in the grand canonical ensemble. While these simulation results confirm 
the existence of the gas-liquid critical point, they fail to provide a conclusive evidence for the nature of 
criticality of the model as well as reliable values for the critical exponents. In addition, the estimates of 
the critical point parameters (.T* = 0.0134 and p* = 0.26) differ from those given in (2), i.e., T* ~ 0.018, 
p* ~ 0.16 (expressed in the same reduced units as in O) which are obtained for smaller system sizes and 
without finite scaling size analysis. In general, the simulation results show a strong dependence of the 
coexistence envelope (its location and shape) on the system size compared to the case of the RPM. The 
gas-liquid phase coexistence in the URPM has been also predicted by the mean-field (MF) like theories 
(71 m, although with the critical point being considerably distant from the simulations. As expected, the 
MF theories predict a classical critical behavior near the critical point. 

Motivated by the above mentioned simulation studies, we address the issue of the gas-liquid criticality 
in the URPM using the theory that exploits the method of collective variables (CVs) (9]|T0). The theory 
enables us to derive, on the microscopic grounds, the effective Hamiltonian and find aU the relevant 
coefficients, including the square-gradient term, within the framework of the same approximation (ill . 
Using this theory, we obtained the effective Hamiltonian of the RPM in the vicinity of the gas-liquid critical 
point that takes into account the contribution from higher order correlations between the positive and 
negative ions (^|T2). The resulting form of this Hamiltonian confirms the fact that its critical behaviour 
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belongs to the universal class of a three-dimensional Ising model. Herein below, we derive the effective 
Hamiltonian of the URPM and, on this basis, analyse the phase and critical behaviour. 

We consider an equimolar mixture of N+ polycations of charge -i-Q and N- (= N+) polyanions of 
charge -Q in a volume V with the total number density p = N/V (N - N+ + N-). The domain of the 
volume V occupied by polyions is denoted by n. Polyions bear an extended charge distribution +Qp(r) 
where the normalized distribution p(r) is supposed to be the same for both species and is given by a 
Gaussian law 


p(r) = 


exp[-r^/{2a^]] 

{2710'^]^^^ 


a - a 12 is the radius of the polyion. The interaction potential between a polyion of species a and one of 
the species jS is as follows O: 


Uapir] 


QaQp 

r 


erf(r/2(T). 


( 1 ) 


This pair potential is finite at a full overlap [erf(r/2CT) ~ r as r ^ 0] ensuring the existence of the thermo¬ 
dynamic limit SO. The Fourier transform of Ua^{r) has the form: 


Uapik) 


47lQaQp 

P 


exp(-P(7^). 


( 2 ) 


The model Q is at equilibrium in the grand canonical ensemble, and hy p- IIk^T we denote the 
inverse temperature (k^ is the Boltzmann constant) and p is the chemical potential. We use the method 
of CVs, with the ideal gas system being a reference system. In this case, one can use the exact functional 
representation of the grand partition function derived for a multicomponent system in fl3l (see equa¬ 
tion (18) in (13l). For the model Q, it can be presented as follows: 


j SiPn J SpQ exp[-^[pN,PQ]), 


(3) 


where the Hamiltonian S^‘[pm,Pq] has the form: 


J^IPn.Pq] 


- J drpiv(r) + ^ J dr [piv(r) -t PQ(r)] In 
^ Jdr [p]v(r)-pQ(r)]ln 


Piv(r)-I-pQ(r) 


H— 
2 


p]v(r)-pQ(r) 


-t ^ Jdrdr' nQQ(|r-r'| 


X PQ(r)pQ(r')-Jdrp]v(r)Viv. 


(4) 


In the above equations, p^ and pq denote the two CVs: p]v(r) = p+(r) -i- p-(r) describing the field of the 
total number density and pQ (r) = p+ (r) - p_ (r) describing the field of the charge density (p+(-) (r), in turn, 
describes the field of the number density of the cations (anions)), uqq is a linear combination of the inter¬ 
action potentials Uap{r]: UQQ{r) - | [u++{r) - 2u+-{r) + u _(r)]. vjv = {v+ + v-)l2, where Valr) = Va-i- 

is the dimensionless chemical potential of the nth species, = Ppa - 3 In Aq., Aq. is the de Broglie thermal 
wavelength and is the self-energy of the nth species = puaair,r]l2 - puol2 with uq - \fna). 
For the model under consideration, UQQ(r) = u++{r) - u _(r) = |u+_(r)| and vjv = v+ = v_ = v. 

Functional integrals (D-d can be given a precise meaning in the case where the domain fl is a 
cube of side L(V = L^) with periodic boundary conditions. This means that we restrict ourselves to the 
fields pA(r) (A = N, Q) which can be written as Fourier series pA(r) = V~^ ZkeA Pk,Aexp(ikr), where A = 
{2jiIL)Z^ is the reciprocal cubic lattice I14II15I . 

Expanding the entropic part of Hamiltonian 0 in powers of p^ and pg (more exactly, in powers 
of deviations of p^ and pg from their MF values), we arrive at the expression similar to that obtained 
in (111 (see equation (7) in (Tl]). The main difference is that in the case of the URPM, the contributions 
from the hard sphere reference system to series expansion coefficients reduce to the ideal gas terms. 
Next, following the procedure developed in (111, one can derive the effective Hamiltonian of the URPM 
near the gas-liquid critical point. Here, we directly use the expressions obtained in fill replacing the 
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contributions from the hard sphere system with the corresponding terms for the ideal gas model. As a 
result, the effective Hamiltonian in the Fourier space representation reads as 




,eff 


ai,oPo,N + 


1 


(^ 2,0 + a2,2) Pk,NP-k, 


N 


+ E —^ 

nHN)’ 


— E ‘^n,0Pki,NPk2,N---Pkn,NSki+...+k„y 

ki,...,kii 


(5) 


where 5kj+.+k„ is the Kronecker symbol. The coefficients in a one-loop approximation have the form: 

1 


^ 1,0 

^ 2,2 


= -Avjv- 


2{N) Es^qW). 




1 


4<iV> 




(6) 

(7) 

(8) 


where Avm-vn- the superscript “id” refers to the ideal gas system. gQ(q) is the screened potential 


gQ(^) = 


PpuQQ(q) 
l + PpuQQlq) 


(9) 


and the superscript ( 2 ) in equation denotes the second-order derivative of gqiq] with respect to the 
wave vector. The terms involving screened potentials arise from integration in over the charge sub¬ 
system (CVs pq). 

Taking into account and one can obtain explicit expressions for coefficients 


aiQ = -v + v 


,id 


2T* 


oc 

-f 

s/n J 




x/n J T* +2Aq x/n 


a„fi = (-l)"(n- 2 )! 


l-(n-l) 


OO 

/ 


(2477)""1 i/TT^-^x^dx 


n^2, 


0.2,2 


-48q 


OO 

x/nT*^ j 


4 2x^ 


[T*+ 24q 


{T*x^e^^ +24q^/nr 

[24q x/^[3 + 7x'^ + 2x*) 


( 10 ) 


( 11 ) 


-T*x^e^ ( 1 -tx^-I- 2 x^)]dx. 


( 12 ) 


In equations 1 10 -i 12 i, the following notations are introduced: T* - k^T!U q - k^T x/nalQ^, q - npa^lG, 
and X = ka. Coefficient ai,a is the excess part of the chemical potential v, and the equation ai,o - 0 yields 
the chemical potential in the random phase approximation (RPA). 

It should be emphasized that although the original Hamiltonian of the URPM [equation dj] does not 
include direct pair attractive interactions of total number density fluctuations, the effective short-range 
attraction does appear in the effective Hamiltonian. This attractive interaction is the result of taking into 
account the charge-charge correlations through integration over the charge subsystem. Therefore, one 
can state that the mechanism of the gas-liquid separation in the URPM is the same as in the RPM. The 
form of the effective Hamiltonian (D-fD suggests the Ising-like critical behaviour of the URPM. Herein 
below, we use explicit expressions (To)-(T^ for the analysis of the gas-liquid phase separation. 

First, we consider the critical point. At the critical point, the system of equations 


a 2 ,o(Pc, Tc) = 0, fl3,o(Pc, rc) = 0 


(13) 


holds. In addition to equations we solve the equation ^ 4,0 = 0. The loci of the above-mentioned 
equations in the (T* ,p*) plane are shown in figure The solid line corresponding to the solutions of the 
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Figure 1. Loci of equations 02,0 = 0 (solid line), 03^0 = 0 (dashed Une), and 04^0 = 0 (dotted line). Tempera¬ 
ture and density are given in dimensionless unites: T* = k^T s/jial[ 2 Q^) and p* = pcr^. 


first equation in is called a spinodal line. The two lines (solid and dashed) intersect at a maximunt of 
the gas-liquid spinodal yielding the coordinates of the critical point. As is seen, coefficient does not 
vanish at the critical point suggesting that this point is a regular critical point. 

Critical temperature and critical density expressed in reduced units are presented in table[^ In order 
to make some contact with the results obtained for the RPM, hereafter we use the same reduced units 
for the density as in l6|, i.e., p* = pa^ where a is the diameter of the polyion. Our choice of the reduced 
temperature, T* - k^Tluo with uq being the maximum strength of the attractive interaction, coincides 
with that of Gllllellil. In particular, one gets uq-2Q^I{ \/na] for the URPM and uo-Q^lcrm the case of 
the RPM where a is the diameter of the polyion in the former case and the diameter of the hard sphere/ion 
in the latter case. With the critical parameters defined in such a way, the critical temperature of the URPM 
appears to be smaller than the critical temperature of the RPM obtained in the same approximation, 
whereas the critical densities of the both models are nearly the same. As expected, our estimates of the 
critical parameters of the URPM coincide with the results obtained in |8) in the RPA. It should be noted 
that the RPA, like other mean-field theories (3, predicts a far too high critical temperature and a far too 
low critical density compared with the available simulation data (l][2][6|. Some possible reasons for such 
a situation have been discussed in |7l. 


Table 1. Reduced gas-liquid critical parameters, coefficients of the effective Hamiltonian and the reduced 
Ginzburg temperature for the ultrasoft restricted primitive model (URPM) and for the restricted prim¬ 
itive model (RPM) in the one-loop approximation. Temperature and density are given in dimensionless 
unites: Tc/ wq and = pc<J^ (see the text for details). The data for the RPM are taken from fill . 


Model 


Pc 



difi 

tc 

URPM 

0.06750 

0.0081 

1.0005 

0.2693 

0.1384 

0.0031 

RPM 

0.08446 

0.0088 

1.0758 

0.2570 

0.1752 

0.0053 


To calculate the coexistence curve, first, equating the right-hand side of (To) to zero, we obtain an ex¬ 
pression for the chemical potential. Then, the Maxwell double-tangent construction is employed. Figurej^ 
shows the coexistence curve of the URPM (a solid line) in the {T*,p*) plane. The coexistence curve of the 
RPM (a dashed line) is presented for comparison. As is seen, the shape of the both curves is very similar. 
The inset shows the both coexistence curves in the corresponding-state variables: t* - {T* - T*)/T* and 
Ap* - (p* - p*)lpl. Comparing the two coexistence envelopes in this case, one can see that the URPM 
has a slightly wider coexistence region than the RPM. 

Finally, using equations and j^, we calculate the reduced Ginzburg temperature fc where to - 
{Tq - Tc)ITc (.Tq is the Ginzburg temperature and Tc is the MF critical temperature). According to the 
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Ginzburg criterion (T3, Tq ntarks a lower bound of the temperature region where a MF description is 
self-consistent. For 1 1\ « tc. Ising critical behavior should be exhibited. One expects a crossover from the 
MF behaviour to the Ising-like critical behaviour when the reduced temperature 1 1\ becomes of the order 
of tQ (16). The reduced Ginzburg temperature can be expressed in terms of coefficients of the Hamiltonian 
as follows (TtI : 


fc ■ 


1 ajo 

3271^ a2,ta\2 


(14) 


where a 2 ,f = da 2 fildt \Taking into account 1 11one can obtain for ^ 2 , 


a2,t - 48r]T^ 


ou 

/ 


dx 


+247] ’ 


(15) 


where the same dimensionless units as in ED are used. 



Figure 2. Gas-liquid coexistence curve of the URPM in the {T* ,p*) representation. The coexistence curve 
of the RPM (a dashed line) is presented for comparison (see the text for details). The inset shows the both 
coexistence curves in the corresponding-state variables. 

Substituting the parameters of the URPM critical point in equations l |TD (for n = 4), (Tp and (T4)-(TD, 
one gets the values for the Hamiltonian coefficients and the reduced Ginzburg temperature in the one- 
loop approximation. The results are presented in table[^ The corresponding values for the RPM obtained 
within the framework of the same approximation in (l3 are presented in the table for comparison. As 
is seen, there is only a little difference between the results obtained for the both models. In particular, 
the models have much the same region of the effective density-density attractions: ^/a^la == 0.5. It is 
essential that the reduced Ginzburg temperature obtained for the URPM is close to that found for the 
RPM, i.e., It is worth noting that the reduced Ginzburg temperature of the RPM is about 

20 times smaller than for the purely nonionic model (a one-component hard-sphere square-well model) 

ED 

In summary, we have derived a microscopic-based effective HamUtonian of the URPM in the vicinity 
of the gas-liquid critical point and have found analytic expressions for all the coefficients in the one-loop 
approximation. As for the RPM, the form of the effective Hamiltonian indicates that the critical behavior 
of the URPM belongs to the universal class of a three-dimensional Ising model. Based on this Hamiltonian, 
we have calculated the reduced Ginzburg temperature fc which, according to the Ginzburg criterion, 
marks the onset of the critical region. We have found that tc for the URPM is of the same order as for the 
RPM. In addition, our results have shown that the shapes of coexistence envelope of the URPM and the 
RPM resemble each other when presented in the similar dimensionless units. Our results suggest that the 
gas-liquid criticality in both models, i.e., the URPM and the RPM, should be of the same nature. 


24601-5 














0. Patsahan 


References 

1. Coslovich D., Hansen J.-P., Kahl G., Soft Matter, 2011, 7,1690; doi 10.1039/C0SM01090A 

2. Coslovich D., Hansen J.-R, Kahl G., J. Chem. Phys., 2011,134, 244514; doi 10.1063/1.3602469 

3. Hansen J.P, McDonald I.R., Theory of simple liquids. Academic Press, 1986. 

4. RueUe D., Statistical Mechanics, Rigorous Results, W.A. Benjamin, New York, 1969. 

5. Fisher M.E., Ruelle D., J. Math. Phys., 1966, 7, 260; doi 10.1063/1.1704928 

6 . Caillol J.-M., Levesque D., J. Chem. Phys., 2014,140, 214505; doi 10.1063/1.4880239 

7. Nikoubashman A., Hansen J.-P, Kahl G., J. Chem. Phys., 2012,137, 094905; doi 10.1063/1.4748378 

8 . Warren P.B., Masters A.]., J. Chem. Phys., 2013,138, 074901; doi 10.1063/1.4791635 

9. Yukhnovskii I.R., Holovko M.F., Statistical Theory of Classical Equilibrium Systems, Naukova Dumka, Kiev, 1980 
(in Russian). 

10. Patsahan O., Mryglod I., Condens. Matter Phys., 2006, 9, 659; doi 10.5488/CMP.9.4.659 

11. Patsahan O.V., Phys. Rev. E, 2013, 88, 022102; doi 10.1103/PhysRevE.88.022102 

12. Patsahan O.V., Mryglod I.M., J. Phys.; Condens. Matter, 2004,16, L235; doi 10.1088/0953-8984/16/16/L01 

13. Patsahan O., Mryglod I., Condens. Matter Phys., 2012,15, 24001; doi 10.3488/CMP.15.24001 

14. Patsahan 0., Mryglod I., Caillol J.-M., J. Phys. Stud., 2007,11,133. 

15. Goldenfeld N., Lectures on Phase Transitions and the Renormalization Group, Addison-Wesley, New-York, 1992. 

16. Chaikin P.M., Lubensky T.C., Principles of Condensed Matter Physics, Cambridge University Press, 1995. 

17. Fisher M.E., Lee B.P, Phys. Rev. Lett., 1996, 77, 3561; doi:10.1103/PhysRevLett.77.3561 


Kpi/iTi/iHHa TOHKa ras-piAi/iHa y/ibTpaM'nKOi o6Me>KeHoT 
npi/iMiTi/iBHOi MOAe/ii: aHa/iiri/mHa reopm 

0. naijaraH 

iHCTi/iTyr (j)i3MKi/i KOHfleHCOBaHMX CMCreM HAH yKpai'Hi/i, ey/i. CBCHqiqbKoro, 1 , 79011 m. /IbBiB, VKpaiHa 

BuBHacTbCti KpuTMHHiCTb tas-pifli/iHa ynbtpaM'tiKOi o6Me>KeHOi npuMin/iBHOi MOfle/ii (URPM) noaieaeKrponiriB, 
Bi/iKopucTOByrotM reopiK), tqo SasyerbCti Ha Merofli KoneKTitBHi/ix smIhhux. fl/ia qiei MOfleni BHBefleno ecfieKTi/iB- 
HHH raMinbTOHiaH i orpMMaHO tiBHi Bnpa3n BCix KoecjjiqiCHTiB B OflHonerneBOMy Ha6nn>KeHHi. Ha ochobI 
qboro raMi/ibTOHiaHy aHaniayeTbCti (fiaaoBa i KpmMHHa noBefliHKa. Orpi/iMaHi peaynbrarM CBi^HaTb npo re, uto 
npwpofla Kpi/iTMHHOcri raa-piflitHa b URPM e raKOio x, ak i b oOMexeniit npi/iMiTi/iBHiit MOfleni, tixa BK/irotae 
TBepflMit Kop. 

KnioHOBi ctioBa: nonieneKTpojiirin, y/ibrpaM'fiKa o6Me>KeHa npi/iMiri/iBHa MOfle/ib, Kpmi^HHa ro^Ka raa-pifii/ina, 
ecjieKT^iBHi^a raMijibTOHiaH 
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